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The damping of spin waves parametrically excited in the magnetic insulator Yttrium Iron Garnet (YIG) 
is controlled by a dc current passed through an adjacent normal-metal film. The experiment is performed 
on a macroscopically sized YIG(100 nm)/Pt(10 nm) bilayer of 4 x 2mm^ lateral dimensions. The spin-wave 
relaxation frequency is determined via the threshold of the parametric instability measured by Brillouin light 
scattering (BLS) spectroscopy. The application of a dc current to the Pt film leads to the formation of a 
spin-polarized electron current normal to the film plane due to the spin Hall effect (SHE). This spin current 
exerts a spin transfer torque (STT) in the YIG film and, thus, changes the spin-wave damping. Depending 
on the polarity of the applied dc current with respect to the magnetization direction, the damping can be 
increased or decreased. The magnitude of its variation is proportional to the applied current. A variation in 
the relaxation frequency of ±7.5% is achieved for an applied dc current density of 5 • 10^^ A/m^. 


The injection of a spin current into a magnetic film 
can generate a spin-transfer torque (STT) that acts on 
the magnetization collinearly to the damping torquei~— 
Thus, it can be used for tuning of the damping of a mag¬ 
netic film^— as well as for the excitation of magnetiza¬ 
tion precession in the filmA”— In the first experimental 
realizations, a dc charge current was sent through an ad¬ 
ditional magnetic layer with a fixed magnetization direc¬ 
tion in order to generate a spin-polarized current 
A different way to generate a spin current is based on 
the spin Hall effect (SHE) ^^d^ caused by spin-dependent 
scattering of electrons in a non-magnetic metal with 
large spin-orbit interaction^^i^iii One of the advantages 
of the SHE is that it does not require a dc current in 
the magnetic layer and, thus, allows for the applica¬ 
tion of a STT to a magnetic dielectric such as yttrium 
iron garnet (YIG)^ which is of particular interest for 
magnon spintronicsi^'^— due to its extremely small damp¬ 
ing parameter Moreover, a great advantage of the 
SHE is that a STT can be applied not only locally but 
to a large area of a magnetic filrr^ and, thus, can be 
potentially used to compensate the damping in a whole 
complex magnonic circuit 

Up to now, the auto-oscillatory regime and the mag¬ 
netization precession generation has only been reached 
in patterned structures Nonlinear multi-magnon 
scattering phenomena are assumed to be the reason 
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that disturbs the generation process in un-patterned 
stuctures. At the same time, nonlinear scattering 
should not affect the SHE-STT-based damping compen¬ 
sation since spin-wave amplitudes in this case are much 
smaller. However, most experimental studies concern¬ 
ing this phenomena using metallic samples^i^ as well 
as YIG structuresiii^ were performed with laterally- 
confined nano- or micro-structures. 

Here, we use a YIG/Pt bilayer of macroscopic size to 
investigate SHE-STT damping variation. The measured 
variation of the damping is proportional to the applied 
dc current and no influence of multi-scattering magnon 
processes is observed. 

The investigated sample consists of a YIG/Pt bilayer 
with macroscopic lateral dimensions of 4 mm x 2 mm, 
and film thicknesses of tyiG = 100 nm and tpt = 10 nm. 
The YIG film is grown by liquid phase epitaxy on a 
gadolinium gallium garnet (GGG) substrate and the Pt 
film is deposited afterwards using plasma cleaning and 
RE sputtering, as described in Ref. Measurements of 
the ferromagnetic resonance (EMR) linewidth and of the 
inverse SHE induced by spin pumping in a wide frequency 
range yield the Gilbert damping parameters ayiG = 
(1.3 ±0.1) • 10“^ for the bare YIG film and GyiG/Pt = 
(5.3 ± O.I) • 10“^ for the YIG/Pt bilayer. In addition, a 
saturation magnetization jaoMs = (173 ± I) mT, an in¬ 
homogeneous broadening /hoAIIq = 0.26 mT, a resistiv¬ 
ity of the Pt film ppt = 1.475 • 10“^ Um, and an effective 
spin mixing conductance = 3.68 • 10^^ m“^ are deter¬ 
mined at room temperature. 

The investigation of the SHE-STT damping variation 
is based on the analysis of the threshold of the parametric 
instability. Eigure [1] shows a sketch of the experimental 
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FIG. 1. (color online) Sketch of the sample and the setup in 
the parallel pumping geometry used for the parametric exci¬ 
tation of spin waves and their detection by BLS spectroscopy. 
The externally applied biasing field and the pumping field 
in the examined area of the sample are both along the x- 
direction, the charge current is applied in-plane along the y- 
direction. 


arrangement, in which the parametric instability thresh¬ 
old is measured in the parallel pumping geometry j^^i^^i^^ 
An external biasing field ii^appi magnetizes the YIG film 
in-plane along the x-axis. The parametric excitation of 
spin waves is achieved by an alternating pumping mag¬ 
netic field hp oriented parallel to i^appi- For this purpose, 
a microwave signal at a fixed frequency of /p = 14 GHz 
is applied to a 50 pm wide Cu microstrip antenna placed 
on top of the sample. A 10 pm thick polyethylene inter¬ 
layer separates the antenna from the sample electrically. 
In our experiment, spin waves at half of the pumping 
frequency (/sw = /p/2 = 7 GHz) are excited as soon 
as the pumping field amplitude hp overcomes a critical 
threshold value hth- The detection of these spin waves 
is realized by means of BLS spectroscopy^^i The incident 
probing laser beam in our experiment, which accesses 
the YIG film through the GGG substrate (see Fig. [1]), 
is always perpendicular to the film plane and only spin 
waves with wavenumbers in a range of k ^ lO^rad/cm 
are detected^^^ Since the threshold /ith depends on the bi¬ 
asing field, both parameters, iLappi and hp, are varied in 
each measurement. Furthermore, An wires are mounted 
to the edges of the sample by silver conductive adhesive 
to apply an in-plane dc current to the Pt film along the 
^-axis. Since the dc current is perpendicular to iLappb 
the SHE generates an out-of-plane spin current (along 
the 2 ;-direction) in the Pt film that exerts a STT on the 
YIG magnetization at the YIG/Pt interface. A maximal 
current of 1 A is used, corresponding to a current density 
of jc = 5 • 10^^ A/m^. In order to reduce the influence of 
Joule heating in the Pt film, the experiment is performed 
in the pulsed regime with a pulse duration of 10 ps and 
a repetition time of 1 ms. All measurements presented 
here are performed at room temperature. 

Figures [2](a)-(c) exemplary show BLS intensities of 
spin waves at /sw = 7 GHz when the pumping field 
hp and the biasing field iJappi (in +x-direction) are 
swept, measured for three different current densities jc 



FIG. 2. (color online) BLS intensity of parametrically 
pumped spin waves at 7 GHz while scanning the biasing field 
Happi and the pumping field hp, for different applied cur¬ 
rent densities (a) jc = +5-10^°A/m^, (b) jc — 0, (c) 
Jc = — 5 • 10^° A/m^. The dashed white line in (b) represents 
the butterfly curve for spin waves with k < 10^ rad/cm. 


(in +^-direction). The applied current densities are (a) 
jc = +5 • 10^^ A/m^, (b) reference measurement jc = 0, 
(c) jc = “5 • 10^^ A/m^. Blue colored areas in the inten¬ 
sity graphs correspond to the case hp < /ith when no spin 
waves are parametrically excited. 

The density n{k,t) of parametrically excited magnons 
of a certain A?-vektor at the time t in the vicinity of the 
threshold is given b y^^i^^ 


n{k,t) = no{k) ■ exp ^-2 (T{k) - |C(fc)|//o^p) t , 


( 1 ) 


where no is the initial magnon density (thermal level), 
F is the relaxation frequency, V{k) is the coupling pa¬ 
rameter between magnons and the pumping field, and /io 
is the vacuum permeability. The magnon density grows 
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FIG. 3. (color online) (a) Obtained Hfmr values (in +x- 
direction) at /ith,min as function of jc. The solid lines are fits 
to Eq.m (b) The induced Oersted field (in +a:-direction) and 
the variation of the saturation magnetization as functions of 
jc, when i^FMR > 0 (in +x-direction). 


exponentially, as soon as the pumping field overcomes a 


. / r(fc) 

m < -= 


critical threshold value of hr^ > /ith = min . 

^ Ooimij 

for one mode with the wavevector k. For a fixed pumping 
frequency, the /ith value strongly depends on the biasing 
field i^appU which determines the wavevector of the avail¬ 
able spin waves. The diagram of /ith vs. iTappi, called 
butterfly curved exhibits a minimum /ith,min at the res¬ 
onant field i^FMR which corresponds to the excitation 
of spin waves with /c ^ 0. The white dashed line in 
Fig. [2Kb) represents only the part of the butterfly curve 
(for the BLS accessible wavenumbers) in a small range 
of the biasing field around its minimum. For the further 
study, only the variations of /ith,min and related 

to k ^ 0 spin waves, are analyzed as functions of the 
applied current density jc. 

First, the variation of the resonant biasing field F^fmr 
corresponding to the minimal threshold pumping field 
^th,min IS investigated. Figure Ha) presents all H^mr 
values obtained from BLS measurements (as shown, e.g., 
in Figs. [ 2 Ka)-(c)) for both directions of the biasing field 
and both directions of the current. iLFMR evidently shifts 
always to values of higher magnitude with increasing jc, 
regardless of the current direction, but the magnitudes 
of the shift are different for opposite current directions. 
The observed behavior can be understood by taking into 
account the contributions of two effects: (i) The dc cur¬ 
rent in the Pt induces an Oersted field iLind collinear to 
^appb which is proportional to jc- Thus, the total mag¬ 


netic field in the YIG film reads 

-f^total — -f^appl T -f^ind — -^^appl H ' jc, (2) 

where the proportionality constant a is introduced as a 
fitting parameter accounting for the sample geometry, 
(ii) Joule heating in the Pt film leads to a temperature 
increase in the bilayer, and decreases the saturation mag¬ 
netization of the YIG film following 

M,(ie) = M,(0) ■ (1 - /? ■ at) = M,(0) ■ (1 - & ■ ic'). (3) 

Here, /io^s(O) = (173 ± 1) mT is the initial saturation 
magnetization at room temperature (300 K) at jc = 0. 
The factor /3 = (2.2 ± 0.1) • 10“^ accounts for the 
change of Mg with a temperature change of AT ^ which 
in turn is assumed to depend on the applied current as 
AT = h/p ' jc^ by introducing the fitting parameter h. 
Using Eq. O Eq. [3] and the Kittel equation {k 0) 

/sw = 7M0\/(77 fMR + 77ind)(77FMR + TJind + ^s), the 

Hymr values in Fig. |3Ka) can be fitted by 

±/^o|77fmr| = -CL • jc T -/^o^s(O) • (1 - 6 • jc^) 

(t)" + ^ }■?))"• * * 

where 7 = 28GHz/T denotes the gyromagnetic ra¬ 
tio. The sign of ±//o|77fmr| indicates the field polar¬ 
ity (” + " corresponds to the +x-direction). The solid 
lines in Fig. |3Ka) represent the fit according to Eq. |4] 
yielding the fitting parameters a = 8.0 • 10 “^^ Vs/A and 
6 = 1.6 • 10 “^^m^/A^. Based on the values of a and b 
the variation of the Oersted field /io^Jind the magne¬ 
tization fioMs with jc are shown in Fig.jSKb) for the case 
of jiioHpMR > 0. The highest applied current density of 
jc = ±5 • 10 ^^ A/m decreases by ~ 7.0 mT, yield¬ 

ing a corresponding temperature increase of AT « 18 K 
for the YIG film. The change of /ioTJind for the maximal 
jc is ~ ±0.4 mT. This magnitude agrees with the the¬ 
oretically expected Oersted field given by (/io/2)jctpt ~ 
±0.3 mT, if the sample width is considered to be much 
larger than the film thickness tyiG- This suggests that 
no STT-based variation of the saturation magnetization 
of YIG is observed in our experiments^ 

Figure HKa) presents the minimal threshold pumping 
field /ith,min as function of jc, extracted from data such 
as shown in Figs. |2Ka) -(c). The shift of /ith,min obviously 
depends on the orientation of jc with respect to iJappi 
and can be understood in terms of the SHE-STT effect. 
According to our experimental geometry, the direction of 
the SHE-generated spin current in the Pt film is given by 
the unity vector {jc x ffappi)/|jc x ffappi|<^ If the spin 
current, for example, is oriented in ± 2 )-direction (from Pt 
to YIG), the damping in the YIG film and consequently 
^th,min dccreascs (compare Fig. l2Kb) with (c)). In the 
case of a reversed spin current direction (from YIG to 
Pt), htii^rnin mcreascs (compare Fig. |2Kb) with (a)). The 
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Current density j, (x10^° A/m^) 


FIG. 4. (color online) (a) Measured minimum threshold val¬ 
ues /ith,min extracted from BLS as a function of jc for positive 
and negative i^FMR. (b) Relative change of the relaxation 
frequency according to Eq. [G] The solid lines are linear fits. 


relaxation frequency F of spin waves with /c ^ 0 as a 
function of jc read o ^ ^ 

TT 

r(4c) = ' Mo-^s(jc) ' M0^th,min(4c): (b) 

^ Jsw 

and the relative variation of F is given by 

^ r(ic) ~ r(Q) _ /?'th,min(ic)-^s(ic) _ /gx 
r(0) " F(0) " /lth,min(0)Ms(0) ' ^ ^ 

Figure [11(b) shows the dependence of the relative vari¬ 
ation of F as a function of jc- These values are ob¬ 
tained according to Eq. [S] using the /ith,min(jc) values 
from Fig. dja), as well as Eq. [3] and the fitting param¬ 
eter b for Ms(jc)- The relative variation of the relax¬ 
ation is proportional to the current within the error bars, 
which is illustrated by the linear fit in Fig |4l(b). F is 
changed by approximately 7.5 % for the highest applied 
current densities of jc = ±5 • 10^^ A/m^ in our experi¬ 
ment. In order to estimate the required critical current 
density for the complete compensation of damping, and 
for the triggering of auto-oscillations in our experiments, 
the extrapolation of the linear fit to r(jfc) = 0 is per¬ 
formed, yielding a value of ^ • lO^^A/m^. This 

value agrees well with the results obtained with patterned 


structuresand with the theoretically estimated value 
of ^ A/m^ calculated on the basis of an¬ 

alytical expressions provided in Ref. 0. For this calcula¬ 
tion, the parameters of our system indicated above, the 
transparency of the YIG/Pt interface^ T = 0.12 (esti¬ 
mated for a spin diffusion length of A = 3.4 nm) and a 
spin Hall angle of ^sh = 0.056 are used^^ Heating ef¬ 
fects are neglected. In addition, reference measurements 
are performed in a geometry with jc parallel to i^appi, so 
that the SHE is eliminated. As a result, no variation of 
F with jc is observed within the error bars (not shown). 

In summary, the threshold of the parametric excita¬ 
tion of spin waves in a macroscopic YIG/Pt bilayer is 
detected by BLS spectroscopy. A change of the damp¬ 
ing in the YIG due to SHE-STT is observed when a dc 
current is passed through the adjacent Pt film. Based on 
our calculations, including current induced Joule heat¬ 
ing and Oersted fields, the damping is found to change 
linearly with the current. The linearity suggests that 
the role of nonlinear multi-magnon scattering processes 
is negligible for the damping variation in the analyzed 
range (below 10%). The maximum used current density 
of jc = ±5 • 10^^ A/m^ results in a damping variation of 
±7.5%. Thus, the complete damping compensation in 
the system is estimated at ^ • 10 ^^ A/m^ which 

is in agreement with the theoretical estimation. 
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